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ABSTRACT 



Introduction: Auditory verbal hallucinations (AVH) in schizophrenia (SZ) have been proposed to result from 
abnormal local, interregional and interhemispheric integration of brain signals in regions involved in language pro- 
duction and perception. This abnormal functional integration may find its base in morphological abnormalities. 
Structurally, AVHs have been frequently linked to abnormal morphology of the superior temporal gyrus (STG), 
but only a few studies investigated the relation of hallucination presence with both whole-brain gray matter 
(GM) and white matter (WM) morphometry. 

Methods: Using a unified voxel-based morphometry-DARTEL approach, we investigated correlates of AVH pres- 
ence in 51 schizophrenia patients (20 non-hallucinating [SZ — ], 31 hallucinating [SZ + ]), and included 51 age 
and sex matched healthy participants. Effects are reported at p < .05 FWE corrected. 

Results: Patients showed lower GM volume of the left STG than controls, irrespective of AVH presence. In addition, 
SZ + showed lower GM volume of the left inferior frontal and right parahippocampal gyrus, and higher WM vol- 
ume of the left postcentral and superior parietal lobule than controls. Finally, volume of the putamen was lower 
in SZ + compared to SZ — . No effects on corpus callosum morphometry were observed. Delusion severity, general 
positive and negative symptomatology illness duration, and medication status could not explain the results. 
Discussion: Results suggest that STG GM abnormalities underlie the general susceptibility to experience psychotic 
symptoms and that additional abnormalities in a network of medial temporal, ventrolateral, putaminal, and parietal 
regions related to verbal memory and speech production may specifically increase the likelihood of experiencing 
AVH. Future studies should clarify the meaning of morphometry abnormalities for functional interregional 
communication. 

© 201 3 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Auditory verbal hallucinations (AVHs) are the predominant type of 
hallucinations experienced by patients with schizophrenia and are 
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characterized by complex auditory phantom perceptions of verbal or 
non-verbal content (Andreasen and Flaum, 1991 ). They have been pro- 
posed to result from abnormal inter-regional and inter-hemispheric in- 
tegration of brain processes (Stephan et al., 2009) in regions associated 
with language and memory processing, including the inferior frontal, 
parahippocampal, and temporo-parietal gyri (Curcic-Blake et al., 2012; 
Dierks et al., 1999; Diederen et al., 2010, 2012; Jardri et al., 2011; 
Vercammen et al., 2010). The abnormal functional integration of 
frontal-temporal language production and perception areas (Curcic- 
Blake et al., 2012) may link to local cortical gray matter (GM) abnormal- 
ities and abnormalities in the pathways that connect these language 
perception and production areas, i.e. the white matter (WM) tracts. 

The superior temporal gyrus (STG) appears to be the most consis- 
tently reported GM abnormality associated with AVH as shown in a 
number of manual tracing studies specifically studying the relation 
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between STG morphometry and AVH (e.g. (Barta et al., 1990; Onitsuka 
et al., 2004; Rajarethinam et al., 2000; Sun et al., 2009)), a finding that 
is supported by two recently performed meta-analyses based on voxel 
based morphometry (VBM) studies related to AVH severity (Modinos 
et al., 2013; Palaniyappan et al., 2012). However, only a limited number 
of independent VBM studies investigated the relations of AVH presence 
and whole-brain morphometry in samples with an adequate size and 
distribution of severity (see (Modinos et al., 2013)). Moreover, only 
one study investigated effects of AVH presence while controlling for 
general disease related factors by including a healthy control group 
(Shapleske et al., 2002). This study however did not find support for the 
specific involvement of STG morphology in AVH. Finally, the majority of 
studies failed to control for other symptoms of reality distortion, includ- 
ing delusions. Therefore, further whole brain studies are needed before 
a definite conclusion on the specific and restricted involvement of the 
STG in the structural pathology underlying hallucinations can be drawn. 

Next to the STG, AVH severity or occurrence has been related to 
morphometric abnormalities of other cortical regions associated with 
language and self-referential processing, including lateral and medial 
temporal-, insular-, inferior frontal-, and medial parietal regions 
(Shapleske et al., 2002; Garcia-Marti et al., 2008; Modinos et al., 2009; 
Nenadic et al., 2010; O'Daly et al., 2007; Shin et al., 2005). Moreover, 
AVHs have been proposed to result from hampered interregional 
communication that finds its base in white matter (WM) volume and 
diffusion abnormalities, thought to reflect the propensity to transfer 
information between regions and hemispheres. Abnormalities in WM 
related to hallucinations have been observed in the temporal, parietal, 
and frontal regions (O'Daly et al., 2007; Shin et al., 2005; Hubl et al., 
2004; Lee et al., 2009; Shergill et al., 2007) and corpus callosum (Hubl 
et al., 2004; Knochel et al., 2012; Makris et al., 2010; Mulert et al., 
2011; Plaze et al., 2011). However, both GM and WM abnormalities 
have been reported with large variability, that could result from variable 
clinical characteristics ( e.g. including only patients with persistent AVHs 
(Garcia-Marti et al., 2008; Modinos et al., 2009; O'Daly et al., 2007) or 
inner vs. outer AVH experience (Plaze et al., 201 1 ) or methodological 
issues such as small (sub)sample size (<20) (Garcia-Marti et al., 2008; 
Mulert et al., 2011; Plaze et al., 2011; Neckelmann et al., 2006)). The 
use of different imaging modalities and analysis procedures (e.g. differ- 
ent segmentation/parcellation strategies), and not taking into account 
total brain volume and severity on other symptom dimensions, includ- 
ing delusional behavior may further contribute to inconsistent results. 

We aimed to test for both GM and WM relations with AVH presence, 
while controlling for delusion severity, total positive symptomatology, 
total negative symptomatology, chronicity of symptoms and type of 
medication used. We investigated whole-brain morphometry in 51 
schizophrenia patients (20 non-hallucinating patients and 31 hallucina- 
tion patients, matched for severity on other symptom dimensions (i.e. 
negative and general symptomatology) and included a sample of 51 
healthy individuals, matched for age and sex. We used a unified voxel- 
based morphometry approach following the Diffeomorphic Anatomical 
Registration Through Exponentiated Lie algebra (DARTEL; (Ashburner, 
2007)), and hypothesized lower GM volume of the STG, inferior frontal 
gyrus (IFG), parahippocampal gyrus, and medial parietal regions specifi- 
cally related to AVH occurrence. With respect to WM, we expected abnor- 
mal morphometiy of the corpus callosum and temporo-parietal region to 
be associated with presence of hallucinations. 

2. Methods 

2.1. Sample characteristics 

Structural MRI data of 102 participants who participated in neuro- 
imaging studies at our center between January 2008 and December 
2010 were pooled for the current analysis (See Table 1 for clinical 
characteristics). Patients were included when complete Positive and 
Negative Symptoms Scale (PANSS; (Kay et al., 1987)) data were 



available from studies conducted at our center (van der Meer et al., 
2013a; Liemburg et al., in preparation; vd Velde et al, in preparation) 
and when data was not included in a previous report on GM correlations 
of AVH severity (Modinos et al., 2009). Studies ran between 2008 and 
2010 at the Neuroimaging Center at the University Medical Center 
Groningen (UMCG) and used the same T1-3D imaging sequence. All 
patients (n = 51) met DSM-1V criteria for schizophrenia. Diagnosis 
was confirmed using the Schedules for Clinical Assessment in Neuro- 
psychiatry or MINl-international neuropsychiatric interview (Sheehan 
et al., 1 998). Forty-one patients were using stable doses of antipsychotic 
medication (classical antipsychotic agents [n = 3]: haloperidol, per- 
phenazine; Atypical antipsychotic agents [n = 38]: risperidone, cloza- 
pine, olanzapine, quetiapine, and/or aripiprazol) and ten patients were 
not using antipsychotic medication at time of inclusion. Healthy con- 
trols (n = 51) were recruited through advertisements and were 
screened to have no Axis-I diagnosis lifetime. Written informed consent 
was obtained from all participants prior to participation to the studies, 
after full explanation of the experimental procedure. All studies obtain- 
ed institutional review board approval (UMCG) and were carried out in 
accordance with the Declaration of Helsinki. None of the current struc- 
tural data has been published as a separate report previously; but sum- 
mary results for GM volumes for part of the sample (N = 33) were 
included in a recent meta-analysis (Modinos et al., 2013). 

2.2. Symptom profile assessment 

Current symptomatology (for the week prior to scanning) was 
assessed using the semi-structured PANSS interview (Kay et al., 1987), 
rated by two trained interviewers. The PANSS has good inter-rater 
reliability and validity (Kay et al., 1988). The interviewers obtained 
a consensus score for each item on all three subscales (positive symp- 
toms (P), negative symptoms (N), and general psychopathology (G)) 
on a seven point scale indicating severity of the symptom in the 
last week (1 = absent;2 = minimal;3 = mild;4 = moderate;5 = 
moderate severe;6 = severe;7 = extreme). For the purpose of study- 
ing volume as a function of hallucination presence, we chose the P3 
(positive subscale, item 3, hallucinatory behavior) score as indicator of 
hallucination severity. Unfortunately, no other measure that reliably 
quantifies intensity, frequency or content of hallucinations was avail- 
able for all patients. Nevertheless, the P3 item of the PANSS has shown 
good convergent validity with the Psychotic Symptom Rating Scale 
(r = .47; (Haddock etal., 1999); and Hamilton Program for Schizophre- 
nia Voices Questionnaire (HPSVQ; r = .66; (Van Lieshout and Goldberg, 
2007)), especially with the factors 'frequency', 'duration', 'disruption', 
and 'beliefs about the origin of the voices' (Kim et al., 2010; Steel et al., 
2007). Because AVHs and delusions could both be considered symp- 
toms of reality distortion and we aimed to specifically test for effects 
of AVH presence, we controlled post-hoc for delusion severity. Item PI 
(positive subscale, item 1, delusions) has shown high correlations 
(r = .88) with the Delusions Severity subscale of the PSYRATS 
(Schneider et al., 2011) and was used to objectify delusion severity. 

2.3. MRI data acquisition 

Imaging data were acquired using a Philips 3-tesla Intera Magnetic 
Resonance system (Best, The Netherlands), equipped with a standard 
SENSE-8 channel head coil. For each subject, anatomical images were 
obtained using a sagittal 3D gradient-echo Tl -weighted sequence 
(TR = 9 ms, TE = 3.5 ms; matrix 256x256; voxel size: lxlxlmm; 
170 slices, duration: 4 min 11 s). 

2.4. Data analysis 

Demographic and clinical data were analyzed with SPSS 16.0 (IBM 
SPSS Statistics, m=nd). Age and sex were entered as nuisance variables 
in each analysis. 
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Table 1 

Clinical characteristics. SZ: Schizophrenia patients; HC: healthy controls; SZ— : non-hallucinating patients; SZ + : hallucinating patients; PANSS: positive and negative syndrome scale; P3: 
PANSS positive subscale, item 3, hallucinatory behavior; PI : PANSS positive subscale, item 1 , delusions. 



Between-group Between group (inl 

HC) 







SZ 


HC 


X2(d) 


F 


P 


SZ- 


SZ + 


X2 d 


F 


P 


N 




51 


51 








20 


31 








Gender (males/female) 


n 


44/7 


37/14 


2.16 




.14 


17/3 


27/4 


3.03 




.22 


Age 


in years; mean (sd) 
range 


34.04 (11.40) 
18-57 


36.14 (10.93) 
18-58 




.90 


.35 


35.00 (9.69) 
20-54 


33.42 (12.50) 
18-57 




.23 


.63 


Education 


in years; mean (sd) 


12.21 (2.87) 


12.64 (2.66) 




.05 


.83 


12.05 (2.59) 


12.32 (3.08) 




.09 


.76 


Handedness (right/left) a 


n 


44/4 


48/3 


.007 




.71 


19/1 


25/3 


.75 




.69 


Gray matter totals 


in mm3; mean (sd) 


759.31 (8.93) 


744.45 (85.32) 




.004 


.95' 


731.23 (69.32) 


777.43 (80.47) 




2.09 


.13 c 


White matter totals 


in mm3; mean (sd) 


503.57 (62.16) 


495.33(62.37) 




.14 


.71" 


479.76 (55.73) 


518.93 (62.05) 




3.39 


.038 c 


csf totals 


in mm3; mean (sd) 


504.95 (62.59) 


497.05 (62.57) 




.11 


.75'- 


481.78 (56.00) 


519.90 (62.88) 




3.18 


.046 c 


Gray + white matter totals 


in mm3; mean (sd) 


1262.88 (112.70) 


1239.79 (133.24) 




.022 


,88 c 


1210.99 (102.72) 


1296.36 (107.39) 




3.15 


.047 c 






















1 x 2 ANOVA 


Age of onset b 


in years; mean (sd) 


25.20 (8.19) 










24.18 (5.750 


25.85 (9.46) 




.43 


.52 


Illness duration b 


in years; mean (sd) 


8.77 (8.71) 










10.41 (7.91) 


7.74 (9.16) 




.98 


33 


PANSS positive symptoms 


mean (sd) 


11.48 (4.19) 










10.84 (4.82) 


11.87 (3.78) 




.71 


.41 


(minus P3) 


range 


7-25 










1-2 


3-5 








PANSS negative symptoms 


mean (sd) 
range 


14.37 (4.18) 
7-24 










12.65 (4.46) 
7-24 


14.84 (3.98) 
9-22 




.99 


.33 


PANSS generalized 


mean (sd) 


29.84 (7.29) 










27.40 (7.16) 


31.42 (7.05) 




3.91 


.054 


psychopathology 


range 


18-49 










18-49 


20-45 








PANSS P3 (hallucinations) 


mean (sd) 
range 


2.55 (1.47) 
1-5 










1.35 (.49) 
1-2 


3.90 (.83) 
3-5 




153.73 


<.001 


PANSS PI (delusions) 


mean (sd) 
range 


2.90 (1.45) 
1-6 










2.25 (1.68) 
1-6 


2.74(1.32) 
1-5 




1.36 


.25 


Classical antipscyhotic use 


n 


3 










1 


2 


2.22 




33 


Atypical antipsychotic use 


n 


38 










17 


21 








No medication 


n 


10 










2 


8 









a Handedness information was not available of 3 patients. 

b RELIABLE information on age of onset and thus illness duration was missing of 7 patients (3 SZ — and 4 SZ+). 
c Corrected for age & gender. 



2.4.1. Preprocessing of imaging data 

Imaging data were analyzed using Voxel Based Morphometry 
(VBM), following the DARTEL approach (Ashburner, 2007) using Statis- 
tical Parametric Mapping software (SPM8; www.fil.ion.ucl.ac.uk/spm) 
implemented in Matlab 7.10.0.499 (The Mathworks inc. Natick, MA). 
VBM-DARTEL preprocessing included: 1 ) manual reorientation of the 
images to the center point of the anterior commissure; 2) segmentation 
of the images into white matter (WM), gray matter (GM), and cerebro- 
spinal fluid (CSF) using the standard unified segmentation option im- 
plemented in SPM8. Image segmentation included an intensity non- 
uniformity correction to account for smooth intensity variations caused 
by gradient distortions (setting: very light regularization); 3) applying 
the DARTEL approach for registration, normalization, and modulation; 

4) further normalization of the GM and WM images to MNI space; and 

5) smoothing of the GM and WM images with an 8 mm full width at 
half maximum Gaussian kernel to increase signal to noise ratio. In the 
resulting images, each voxel represented an absolute amount of brain 
volume, equivalent to the brain volume per unit prior to normalization. 
To achieve maximal sensitivity, to optimize voxel residual smoothness 
estimation and to exclude false positives in non-GM (or non-WM) tis- 
sue, voxel-wise comparisons were masked using a comparison- 
specific explicit optimal threshold GM (WM) mask created using the 
Masking toolbox (//wwwO.cs.ucl.ac.uk/staff/g.ridgway/masking/). 

2.4.2. Statistical modeling 

Data were analyzed in the context of the General Linear Model for 
GM and WM images separately. We performed an ANCOVA to test for 
the main effect of group with group (3; healthy controls, no hallucina- 
tions [SZ — ; P3 < 3], hallucinations [SZ +; P3 > 3]) as random factor, 
the GM (WM) density maps as dependent factor, and age and sex as co- 
variates. Total parenchyma (GM + WM totals; globals) were entered 
by means of proportional scaling. Patients with a P3 score of 1 or 2 
were classified as having no hallucinations (SZ — ), whereas patients 



with a P3 score of 3 and higher were classified as 'hallucinators' (SZ + ). 
To test for the possible confounding influence of other clinical variables, 
we extracted volumes of clusters showing a main effect of group at 
p < .005 and exported these to SPSS. In a first post-hoc linear regression 
model (method enter), we tested whether PI (delusions) and positive 
symptoms (minus P3 & PI) and total negative symptomatology ex- 
plained any additional variance or could explain the effects instead. In 
a second model, we entered years since onset of first psychotic episode 
(i.e. chronicity) to the model. In a third model, type of anti-psychotic 
used were stepwise entered into the model by means of 3 dummy vari- 
ables coding for either use of 2nd generation antipsychotic use (n = 38), 
classic anti-psychotics (n = 3), or no medication use (n = 10). Unfor- 
tunately, accurate information on dose and duration of medication use 
was not available of all patients, excluding the possibility to control 
for medication-status using dose-equivalences (e.g. chlorpromazine- 
equivalences or the olanzapine-equivalences). Hallucinations yes/no, 
age, sex, and GM + WM totals were first entered to all the models. 

2.4.3. Statistical thresholding 

Significance for demographical and clinical data was set at p < .05, 
two-tailed. 

Main effect of Group (HC, SZ— , SZ+) were explored using an F-tests 
atp < .005 uncorrected. Post-hoc T-tests had to meetp < .05 Family Wise 
Error (FWE) corrected at the cluster level to be considered significant 
using a voxel-wise threshold of p < .005 uncorrected. Non-stationarity 
correction was applied to account for the non-uniformity of the 
smoothness of VBM data (http://fmi"i.wfubmc.edu/cms/software#NS). 

Multiple comparison correction of effects occurring in our a priori re- 
gions of interest was restricted to the spatial extent of a composite ROI 
covering these regions and was constructed based on imaging reports 
that relate to hallucinations as outlined in the introduction. The composite 
ROI encompassed the bilateral IFG, bilateral temporal parietal areas, the 
posterior cingulate cortex (PCC), and the parahippocampal regions. The 
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following labels of the Anatomic Automatic Labeling system implemented 
in the WFU Pickatlas (http://www.nitrc.org/projects/wfu_pickatlas/) 
were used: Fronta Lin f_Oper_L, Fronta l_In f_Tri_L, Fronta l_In f_Oper_R, 
Frontal_Inf_Tri_R, Tempora l_Sup_L, Parietal_Inf_L, Temporal_Sup_R, and 
Parietal_lnf_R, Heschl's_gyrus_L, Heschl's_gyrus_R, Cingulum_Post_L, 
Cingulum_Post_R, ParahippocampaLL, and ParahippocampaLR labels. 
As defining ROIs for WM regions is less straightforward, we used 
bilateral lobar labels covering the whole lobule implemented in the 
WFU pickatlas to define the correction area for white matter effects 
(i.e. bilateral whole temporal lobe for temporal effects, bilateral whole 
parietal lobe for parietal effects, bilateral whole frontal lobe for frontal 
effects and bilateral whole limbic lobe for limbic effects including the 
posterior cingulate cortex). Effects occurring outside these ROIs had 
to meet cluster wise p < .05 FWE whole-brain correction (following 
non-stationarity correction). 

3. Results 

3.1. Sample descriptive 

Patients did not differ from HC in age, years of education, handedness, 
sex, and total GM, total WM, or total GM + WM volume (Table 1). 
Hallucinating patients (SZ + ) did not differ from non-hallucinating 
patients (SZ — ) in age, years of education, handedness, sex distribution, 
age of onset of the disorder, total positive symptomatology (minus P3), 
total negative symptomatology, total general psychopathology, delusion 
severity, and medication status (Table 1). However, SZ+ had higher 
WM and WM + GM matter volume than SZ - (post-hoc t-test: WM: 
p = .035; GM + WM: p = .042; all Bonferroni corrected), but not 
compared to HC. 

3.2. VBM results: gray matter 

A main effect of group was observed in the left STG extending into 
the inferior parietal lobule (1PL), the left inferior frontal gyrus (1FG), 
the bilateral parahippocampal gyrus, bilateral putamen, the medial 
frontal gyrus, the precentral gyrus, and the right lingual gyrus (See the 
Appendix lc, Table A-I for a complete overview and details of the main 
effects). Plotting of parameter estimates guided the post-hoc t-tests 
in these regions and showed that patients were characterized by 
lower volume of the STG, an effect that was independent of hallucina- 
tion presence (See Fig. la). 

Compared to HC, SZ + showed lower volume of the left IFG and right 
parahippocampal gyrus (Fig. lb). The left parahippocampal gyrus was 
observed just subthreshold in this comparison (k = 836, Z = 3.93, 
MM coordinate: [z = -21y = 3z = -29], Pfwe_roi = -085). No 
significant volumetric differences in these regions were observed 
between SZ+ and SZ— , and between SZ — and HC. 

SZ + showed lower volume of the left putamen than SZ — , an 
effect that was observed just sub-threshold in the right putamen 
(k = 1765, Z = 3.89, MNl-coordinate: [z = 29 y = 5 z = - 11], 
PFWE_whoieBrain = -19). As Fig. lcl shows, an inverted u-shape is sug- 
gested in the relation of putaminal volume and hallucination presence, 
as SZ — appeared to show higher putaminal volume than HC. However, 
this difference was not significant. Finally within our composite ROI, 
SZ— showed lower volume of the bilateral PCC compared to SZ+, 
although this effect did not reach significance (k = 591, Z = 3.35, 
MN1 coordinates: [x = -8 y = -36 z = 31] / [x = 9 y = -40 
z = 28], Pfwejoi = -27). 

Results are listed in Table 2. 

3.2.1. Effects of confounders 

Volumes of the significant clusters per participant were exported to 
SPSS to test for confounders. Within patients, volume of the STG was 
unrelated to P3 score (continuous), both across SZ+ and SZ— and 
within SZ + (SZ: (3 = 0.027, p = .78; SZ + : [i = 0.13, p = .35). Also, 



STG volume was unrelated to Positive symptomatology (minus PI & 
P3), Negative symptomatology, years since disease onset, and type of 
medication (all p < .14, p > .28). 

Volume of the IFG and parahippocampal gyrus within SZ + was not 
further correlated with hallucination severity, and volume was unrelat- 
ed to PI, Positive psychopathology (minus PI & P3), Negative symp- 
tomatology, age of onset and medication used (all [i < |.45|, p > .16). 
Finally, the difference between SZ + and SZ — in putaminal volume 
was unaffected when clinical and medication confounders were step- 
wise added to the model. 

3.3. VBM results: white matter 

A main effect of group was observed on WM of the left superior and 
inferior parietal lobule, the superior temporal lobule, and the right mid- 
dle frontal lobule (See the appendix, Table A-l for a complete overview 
and details of the main effects). 

Post-hoc t-test showed that SZ+ had higher WM than SZ- in the 
postcentral gyrus extending into the superior parietal lobule and medial 
precuneus (See Table 2 and Fig. 2a) and subthreshold in the 1PL 
(k = 447, Z = 3.82, MN1 coordinates: [x = -29 y = -46 z = 40], 
Pfwe_roi = .11) and PCC (k = 209/66.Z = 3.42/3.19, MM coordinates: 
[x = -29 y = -46 z = 40], PfWE.Ro/ = -21). The relations of WM 
volume in the superior and inferior parietal lobule were characterized 
by (inverted) U-shapes (See Fig. 2), as SZ — appeared to have lower 
WM of the superior parietal lobule than HC (k = 65, Z = 3.56, MN1- 
coordinates: [x = — 26 y = —40 z = 60]) but higher volume of the 
IPLthanHC(k = 183, Z = 3.33, MN1 coordinates: [x = -48y = -46 
z = 39]. However, these comparisons did not reach significance 
(PiWE.Ro/ = -45 and Pfwe_ro/ = -27, respectively). 

Compared to HC, SZ + had higher WM of the superior temporal lobe, 
adjacent to the GM cluster that was observed in SZ — and SZ + com- 
pared with HC (See Fig. 2b). However, the WM effect did not reach sig- 
nificance (k = 307, Z = 4.07, MNI coordinates: [x = -60 y = -27 
z = 3], Pfwe_ro/ = -23). 

Finally, no effects of hallucination presence were observed in the 
corpus callosum, also not at p < .005 uncorrected. 

3.4. Post-hoc sensitivity analysis 

Patients that demonstrated questionable hallucinatory symptom- 
atology scored a '2' on the PANSS-P3 item. Because of the uncertainty 
in hallucinating behavior, these patients may obscure the relation 
between regional brain volume and absence/presence of AVHs. There- 
fore, we post-hoc repeated the analysis on the volumes associated 
with AVH-presence after excluding seven 2-scoring patients using 
SPSS. However, excluding these seven patients did not influence the 
observations. If any, it made the association with putaminal GM- 
volume and SPL WM-volume stronger. Details on these analyses and 
results can be found in Appendix-lb. 

4. Discussion 

In this study, we investigated GM and WM morphometry related to 
AVH occurrence in 51 patients with schizophrenia (20 without halluci- 
nations [SZ — ]; 31 with hallucinations [SZ + ]) and 51 healthy controls. 
Results indicate that lower volume of the left superior temporal gyrus 
(STG), including the primary auditory cortex is a robust abnormality 
in patients with schizophrenia, which seems unrelated to hallucination 
presence or severity (confirmed by whole-brain continuous linear 
regression analyses; see Appendix-1 c and Fig. A-l ). Abnormal low vol- 
umes of the left inferior frontal gyrus (IFG) and parahippocampal 
gyrus were however specifically related to the presence of hallucina- 
tions. Our results further indicate that the putamen plays a role in the 
structural pathology underlying hallucinations, where low volume is 
characteristic of SZ+ as compared with SZ— . Finally, higher absolute 
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Fig. 1. Gray matter volumetric effects, a) Patients (both hallucinating [SZ + ] and non-hallucinating patients [SZ— ]) showed lower volume of the left superior temporal gyrus (STG) than 
healthy controls [HC], extending into the inferior parietal lobule, no difference in STG volume was observed between SZ + and SZ— (see plotAl). b) SZ+ showed lower volume than HCof 
the parahippocampal gyri (of which only the right parahippocampal effect was significant; see plot Bl) and left inferior frontal gyrus (IFG: see plotB2). Although volumes appeared lower 
in SZ— compared to HC, results were non-significant, c) SZ+ showed lower volume of the bilateral putamen than SZ— , although only the left side survived whole brain correction for 
multiple comparisons. All effects are displayed at p < .005 uncorrected. Plots Al, Bl, B2, and CI represent the mean GM volume of the clusters (corrected for total brain volume, age, 
and gender) revealed by the F-test testing for the main effect of group. The y-axis denotes mean volume in ml. 



whole-brain WM and relative WM volume of the left postcentral and 
parietal lobule was specifically associated with occurrence of hallucina- 
tions. All effects were unrelated to delusion severity, total positive 
symptomatology (minus hallucinations and delusions), negative symp- 
tomatology, general psychopathology, age of onset of the disorder or 
medication status. No effect on corpus callosum morphometry was 



observed. Together our results suggest that morphological abnormali- 
ties in a structural network encompassing the IFG, parahippocampal 
gyrus and putamen in addition to the disease general STG abnormalities 
may underlie the experience of hallucinations. 

The medial temporal gyrus, inferior frontal gyrus and putamen 
abnormalities specific to hallucinations, together with the disease- 
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Table 2 

between group comparisons of regional gray and white matter volume. L/R: left/right hemisphere; BA: Brodmann area; k: cluster size at p < .005, uncorrected, MNI: Montreal Neurological 
Institute; p (whole brain): whole brain FWE corrected p-value at the cluster level after applying non-stationarity correction; p(ROI): FWE corrected p-value at the cluster level for the ex- 
tent of the composite ROI after applying non-stationarity correction. 

Gray matter between-group comparisons 
SZ<HC 



MNI coordinate 



L/R 


Region 


BA 


k (whole brain) 


k (ROI) 


X 


y 


z 


Z-score 


T-value 


p (whole brain ) 


P (ROI) 


L 


Superior temporal gyrus 


41 


2194 


883 


-50 


-28 


15 


3.78 


3.94 


.036 


.043 




Transverse temporal gyrus 


42 






-60 


-12 


10 


3.31 


3.42 








Posterior insula 


13 






-44 


-9 


18 


2.64 


2.70 







SZ+ < sz 



L/R 



Region 



BA 



k (whole brain) 



MNI coordinate 



Z-score 



T-value 



p (whole brain) 



p (ROI) 



Putamen 
Putamen 
Putamen 



2146 



-33 
-29 
-24 



-2 
-9 
-11 



3.47 
3.47 
3.44 



3.59 
3.59 
3.56 



.048 



SZ+ < HC 



MNI coordinate 



L/R 


Region 


BA 


k (whole brain) 


k (ROI) 


X 


y 


z 


Z-score 


T-value 


p (whole brain) 


P (ROI) 


R 


Uncus/parahippocampal gyrus 


28/35 


1427 


707 


21 


6 


-29 


4.55 


4.81 


.12 


.044 


L 


Inferior frontal gyrus 


47 


1369 


548 


-42 


35 


0 


3.86 


4.02 


.06 


.019 




Insula 


13 






-41 


11 


7 


3.34 


3.45 








Inferior frontal gyrus 


47 






-41 


24 


0 


3.20 


3.29 







White matter between-group comparisons 
SZ- < SZ+ 



MNI coordinate 



L/R 


Region 


k (whole brain) 




X 


y 


z 


Z-score 


T-value 


p (whole brain) 


P (.ROI) 


L 


postcentral lobule 


518 


492 


-30 


-46 


61 


3.91 


4.09 


.10 


.029 




superior parietal lobule 






-18 


-51 


61 


3.58 


3.72 








medial parietal lobule/precuneus 






-9 


-58 


52 


3.06 


3.14 







generic STG abnormality, strikingly resembles the network of brain 
regions recently described in meta-analyses of functional studies on 
voice hearing in both psychotic and non-psychotic people (Jardri 
et al., 2011; van et al., 2013). Also, this set of morphological abnormal- 
ities supports the corticostriatal model proposed by Hoffman et al. 
(Hoffman et al., 201 la; Hoffman and McGlashan, 2001 ). In this model, 
increased putamen-IFG connectivity is proposed to be specific for 
SZ+, whereas increased STG-putamen connectivity is common for SZ 
patients irrespective of hallucination presence. We speculate that this 
hyperconnectivity may result from lower cortical volume owing to a 
lower propensity to exert inhibitory control in the communication 
between regions. Lower cortical gray matter may represent a lower 
glial count, lower eel size, or low density of neurons in the putamen, 
both excitatory and inhibitory. Lower glial and neuronal density may 
represent a loss of inhibitory activity, causing disinhibition from one 
region to another and higher connectivity between regions. Along this 
line of reasoning, the lower 1FG and putaminal volume in SZ + patients, 
may result in a lack of down-regulation of the 1FG over the putamen, 
leading to a more pronounced hyperconnectivity between these 
regions in hallucinating patients (SZ + ). However, one might also 
argue that higher volume leads to a higher propensity to communicate 
and thus to higher connectivity, indicating a need for future studies to 
focus on structure-function relations. Our current results reflect the 
commonness of STG abnormalities for SZ patients and specificity of 
putaminal and IFG abnormalities for hallucinating patients. Moreover, 
we found abnormal parahippocampal morphometry to be more abnor- 
mal in hallucinating than in non-hallucinating patients as compared 
with controls. The parahippocampal gyrus has, within the context of 
voice hearing, been linked to memory processing and retrieving content 



from memory (Diederen et al., 2010) and abnormality in these process- 
es may corrupt narratives in language (Hoffman et al., 201 lb). Howev- 
er, this hippocampal region has also been linked to generation of 
conscious preceptions (Behrendt, 2013), binding of information from 
different sensory association cortices (Behrendt, 2013), behavioral 
inhibition and emotion regulation (Sjoerds et al., 2013; Gray and 
McNaughton, 2003), and attentional regulation over sensory inputs 
(Behrendt, 2010), and thereby may play an important role in the 
occurrence and content of hallucinations. Moreover, functional connec- 
tivity patterns of the hippocampus was shown to differentiate those 
experiencing auditory from those experiencing auditory and visual 
hallucinations (Amad et al., 2013), further suggesting the importance 
of the hippocampus for integrating signals from multiple sensory 
cortices for conscious experiencing (Behrendt, 2010). 

Finally, the importance of this regions for experiencing hallucina- 
tions, for altered consciousness, and abnormal integration of informa- 
tion processed in cortical areas for accurate sensory representations, 
and finally note that differential connectivity of this region may link to 
the content of hallucinations (i.e. with or without visual hallucinations 
present). The left IFG has been linked to both voice hearing and voice 
hearing detection in functional neuroimaging studies (Diederen et al., 
2012; van et al., 2013). Our results support the role of the left IFG in 
the primary pathology of hallucinations. Next to its undisputable role 
in language production, the IFG has been linked to attentional and in- 
hibitory control and the selection of goal-appropriate behavior, underly- 
ing adequate emotion and behavioral regulation (Ochsner et al., 2012; 
van der Meer et al., 2013b). Finally, the putamen has been linked to 
the initiation and execution of speech (Price, 2010; Argyropoulos et al., 
2013), integration of auditory and visual information (von Saldern and 
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Fig. 2. White matter volumetric effects, a) Patients (both hallucinating [SZ + ]) showed higher postcentral and superior parietal lobule white matter (WM), extending into the medial 
parietal lobule and inferior parietal lobule, as compared to non-hallucinating patients [SZ— ] (see plot Al). b) SZ+ showed higher superior temporal lobule as compared to HC, adjacent 
to the gray matter (GM) cluster that was lower in volume in patients compared to HC, irrespective of hallucination presence. The GM cluster is depicted in yellow, the WM cluster is 
depicted in blue. All effects are displayed atp < .005 uncorrected. Plot Al represents the mean WM volume of the postcentral/superior parietal lobule cluster (corrected for total brain 
volume, age, and gender) revealed by the F-test testing for the main effect of group. The y-axis denotes mean volume in ml. (For interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this article.) 



Noppeney, 2013), and with the conscious perception of auditory stimuli 
(Mhuircheartaigh et al., 2010), underlining the significance of the puta- 
men for speech perception, production and multimodal integration that 
all have importance for experiencing hallucinations. 

As GM volume of the STG was unrelated to current hallucination 
presence and severity, positive symptomatology in general, and dura- 
tion of illness, we suggest that abnormalities in certain portions of 
the STG may not be truly specific to AVHs, but rather related to the 
schizophrenia disease process in general. The amount of studies that 
related schizophrenia to STG pathology is abundant (Sun et al., 2009; 
Honea et al., 2005) and places the STG at a prominent position within 
the structural pathology underlying the occurrence of psychotic symp- 
toms, although its specific relevance for hallucinations has been report- 
ed as well (Palaniyappan et al., 2012; Neckelmann et al., 2006; Gaser 
et al., 2004). Discrepancies in results could be explained by differences 
in sample size (Neckelmann et al., 2006), number of hallucinating 
patients (Gaser et al., 2004), or a non-optimal distribution of hallucina- 
tion severity to study correlations in a linear fashion (Nenadic et al., 
2010). It is important to note that we do not refute the importance of 
the STG for hallucinations, but propose that the likelihood for 
experiencing hallucinations is increased by the specific involvement of 
a defective network of medial temporal, putaminal and lateral prefron- 
tal regions in addition to general disease-related STG pathology. 

An intriguing open issue is how structural variation influences func- 
tional integration. In attempts to link function to structure it became ap- 
parent that relations are not always mapped one to one: functional 
relations may exist in the absence of structural connections whereas a 
fair amount of information is transferred via indirect structural 



connections (Honey et al., 2009). In the specific case of schizophrenia 
patients, researchers have noted that function-structure relations are 
stronger in patients (van den Heuvel et al., 2013), which would suggest 
that our volumetric results predict decreased connectivity between 
these regions, in line with previous work of our group (Curcic-Blake 
et al., 2012; Vercammen et al., 2010). However, structure-function rela- 
tions have only been studied using diffusion tensor imaging data. So far, 
the predictive value of GM morphometry for functional connectivity 
has, to our knowledge, not been studied in schizophrenia, although in- 
teresting relations have been shown in major depressive disorder (van 
Tol et al., 2013). We here hypothesize that decreased volume of the 
IFG and parahippocampal gyrus and putamen in SZ + together with a 
disease-generic abnormality in the STG hampers the direct connectivity 
between these regions, leading to a disintegration of signal and eventu- 
ally to the experience of hallucinations. 

Related to WM abnormalities, we found no associations with corpus 
callosum volume, but relations between AVH and total WM volume and 
relative WM volume of the postcentral and superior parietal lobule 
were observed, and partly included the arcuate fasciculus. This WM 
tract has been repeatedly linked to AVH in schizophrenia in diffusion 
tensor imaging studies. Our result of increased WM volume fits the 
previous observations of higher fractional anisotropy (FA) (Hubl et al., 
2004; Abdul-Rahman et al., 2012) and higher volume of the arcuate 
fasciculus (Hubl et al., 2010). As higher FA and volume most likely 
reflects a higher propensity to transmit signals between regions 
(Witelson, 1995), increased WM in the region connecting this post- 
central gyrus and medial parietal regions may increase the likelihood 
that information is transmitted between temporo-parietal cortical 
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regions associated with voice hearing (Diederen et al., 2012; Jardri et al., 
2011; Kuhn and Gallinat, 201 0), the generation of inner speech (Shergill 
et al., 2003), and self-referential processing (e.g. the medial precuneus/ 
PCC (Whitfield-Gabrieli and Ford, 2012). However, conflicting results 
have been reported as well. Decreased arcuate fasciculus FA in patients 
with a history of hallucinations as compared to patients without such 
history (Catani et al., 2011), compared to healthy controls (de Weijer 
et al., 2011), or compared non-psychotic voice-hearers (de Weijer 
et al., 2013) has been reported as well. However, these results may 
not specifically apply to the posterior parts of the AF (Catani et al., 
2011), and specificity of effects for hallucinating vs. non-hallucinating 
patients were not explicitly studies (de Weijer et al., 201 1, 2013). More- 
over, we could not support lower FA in the posterior parts of the AF in a 
subsample of the current sample ( Curcic-Blake et al., 201 3 ) . The absence 
of a relation between hallucinations and volume of the corpus callosum 
does not support a central role of this structure in the structural pathol- 
ogy underlying AVH. This is consistent with earlier studies (Rossell et al., 
2001 ; Rotarska-Jagiela et al., 2008; Serpa et al., 2012), although abnor- 
mal diffusion abnormalities have been observed related to hallucina- 
tions (Hubl et al., 2004; Knochel et al., 2012; Makris et al., 2010; 
Mulert et al., 2011). However, current macroscopic observations can 
not fully exclude the possibility that microscopic abnormalities in inter- 
hemispheric communication through the corpus callosum underlies 
hallucinations. 

Some limitations should be mentioned. Although medication use did 
not differ between hallucinating and non-hallucinating patients and we 
investigated if use (yes/no) affected the observed effects, a possible 
effect of dose and duration of usage of antipsychotic medication on 
our results cannot be ruled out (Ho et al., 2011). Unfortunately, dose- 
equivalences of for example chlorpromazine could not be calculated 
because reliable dose information was not available of all patients. 
Moreover, we employed only one measure of hallucination severity; 
thereby we had limited information on frequency, content, duration, 
and form of the hallucinations. Although the majority of patients expe- 
rienced auditory-verbal hallucinations as expressed during the inter- 
view, use of the PANSS gives no specific information on content of 
hallucinations and therefore current results may not be fully specific 
for AVH, but relate to the general presence of hallucinations. Moreover, 
next to reflecting current propensity to hallucinate (last week), the 
PANSS P3 item may additionally reflect efficacy of treatment or levels 
of experienced stress in the past week. Caution should therefore be 
taken when interpreting our results as reflecting trait-or state hallucina- 
tions. Nevertheless, in various studies the P3 item has been shown to be 
a valid measure of AVH presence (Kim et al., 2010; Steel et al., 
2007), and to be highly correlated with more extensive assessments, 
such as the Psychotic Symptoms Rating Scale (PSYRATS; (Haddock 
et al., 1999) and with the factor AVH hallucination duration of the 
Schizophrenia Voices Questionnaire (HPSVQ; (Kim et al., 2010; Steel 
et al., 2007). 

4.1. Conclusions 

In sum, results from this study suggest that structural abnormalities 
in a cortical network including the medial temporal lobe, inferior frontal 
gyrus and putamen in addition to disease-common abnormal superior 
temporal morphology increases the likelihood for experiencing halluci- 
nations. Abnormal WM morphometry of the postcentral and parietal 
lobule may further promote interregional communication subserving 
voice hearing and self-referential processing. Future studies should 
focus on the meaning of CM and WM morphological abnormalities for 
integration of functional networks associated with AVH. 
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